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3UMMARY 


Th.  standard  algorithm  of  th.  revised  slmplax  math,  l-*ith  which 
.  reader  1.  assumed  to  b.  familiar^..  been  designed  co  provide 
th.  Slmpiaat  possible  s.t  of  rules  for  solving  a  ltar  progrMmlng 

pro  ..  in  the  general  case  where  nothing  1.  too*,  or  assumed  about 
the  ran*  of  th.  system.  Bxpen.no.  indicate,  that  this  simplicity 

IhlTTT"  *’  Sal"*d  “  th°  several  iteration, 

-hi  h  tend  to  replace  certain  vectors  over  tod  over.  It  1.  believed 

that  this  tendency  can  be  reduced  by  a  mora  critical  „.iy,ia  of  th. 

objectives  of  th.  procedure  on  each  Iteration.  While  this  compll- 

rU1*“  °f  somewhat, .  the  latter  le  not  a  ..nous 

obstacle  for  large,  high-speed  computers. 

Previous  paper.  Involving  th.  us.  of  a  dual  algorithm  have 
suggested  that  much  Information  1.  inherent  In  a  problem  which  1. 
not  normally  taken  sd  van  tag.  of.  Th,  prasent  proposal  u...  a  nom¬ 
ination  of  th.  normal  tod  dual  algorithms,  with  some  modification.-, 

kMP8  “Par*te  th*  th"'  ««n‘i«ly  distinct  Idea,  with  which 
the  Simplex  method  1.  concerned,  namely,  non-elngularlty  of  a  basis 

matrix,  feasibility  of  th.  solution,  and  optimality  of  th.  solution. 

is  believed  that  to  optimal  solution  will  be  obtained  In  fewer 
iterations  by  this  method. 
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A  COMPOS ITS  SIMPLEX  ALQQRXTHH  II 
Wm.  Orchard— Hays 
I— NOTATION  AND  BACKGROUND. 

*  linear  programing  problem  1.  ooncemed  with  mlnlml.lng  (or 
maximizing,  which  merely  involve.  chang.  of  eigne,  a  lln.ar  fonn 


(1) 


ill*' 


oj  XJ  •rain 


subject  to  the  conditions  that 

(2)  7  a.  .  x.  ■  b.  <  a  1  o 

A  A •  |in 


A  •«  xj  ■  bi  • 


(3) 


X.  >  0  ,  J-1P  n 

J  9  v  ±  9  ^  p  •  •  •  pYl  p 


where  we  assume  (without  loss  of  generality)  that  all  b,  ;>  o. 

In  the  revised  simplex  method  (and  m  particular  here,  the  fom 
of  the  method  using  the  product  fonn  for  the  inverse  of  the  basis) 
the  statement,  (l),  (2)  are  combined  in  a  single  matrix  equation 
formed  of  the  following  columns  (denoted  by  braces): 

Po  "  •••»  0 1 

PJ  "  {*oy  alj'  •••»  amjj  '  J  ■  1,2, ... ,n 
Q  "  (°'  bl*  •••*  bm} 

X  "  {x o'  xl»  •••*  *n} 

which  are  combined  to  form  the  matrix  p  -  Tp  p  p -j  .  .. 

L o'  i>  •••»  *VJ  »nd  the 

equation 


(4)  PX  -  Q 


*  n 
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whence  xc  is  to  be  maximized  subject  to  (3)  and  (4). 

**•  simplex  method  works  solely  with  basic  solutions.  If  p 
does  not  contain  the  (m+l)-order  Identity  matrix  I,  then  a  basic 

~ g-blC  Bolutlon  18  not  m  general  available  to  start  the  Iterative 
procedure  for  maximizing  xQ.  That  Is,  what  Is  desired  Is  the  ex¬ 
plicit  knowledge  of  the  following  quantities: 

(5a)  a  basis  for  (m+l)-space  chosen  from  the  P  and  denoted 
by 

®  m  Q* 4  #  ^4  »  ....  P.*l  .  p  m  P 

Jo  Ji  4  4  0  ' 

(5b)  its  associated  solution  vector  V  -  lv  v  „  7 

l  o'  vi»  vmJ 

such  that  BV  -  Q,  v±  ^  0  for  1  >  0  ;  and 

(5c)  the  inverse  of  B  whose  elements  we  denote  by  0lk,  the 

row.  by  px  -  (plo,  Pa . Pij  ,  i  .  0<1 . . 

«d  th.  columns,  when  necessary,  by  C*  .{pok,  . . 

^  m  ®*i#»s.#m  ;  that  Is,  briefly, 

B  1  -  (Plk)  -  {pQ,  pmJ 

where  the  last  braces  signify  a  column  of  row. 

Besides  the  lack  of  knowledge  of  the  quantities  (5),  it  is  in 
general  not  even  known  whether  or  not  the  matrix  P  has  rank  m+1, 
that  is,  whether  or  not  it  is  possible  to  choose  a  basis  for  (m+l)- 
space  from  the  columns  of  P.  At  the  same  time,  if  certain  of  the 
equations  (2)  were  redundant  (i.e.,  rank  of  P  less  than  m+l),  this 
would  not  preclude  the  possibility  of  a  solution  to  the  problem. 
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perhaps  a  unique  one.  Furthermore,  even  If  the  quantities  (s) 

*"  kno"n'  «  —m  b.  t0  Btmrt  th#  (5) 

the  simplex  method  using  the  product  form  ox  B”"1  m  a  ^ 

“r~- 

:z:: — 

olutlon  or  show  that  non.  existed.  It  i.  assumed  that  the  reader 

:::r:rth  thi-  p~~du~  -  -  - -  -  -  — t  • 

B_1  -  *k  V.  •••  *! 

Wh"r*  th*  S*  ar®  el*msntary  oolumn  matrices.  Cl, 3,5,6] 

ix-disadvantaoes  of  phabb  one  OF  mb  sxkpxjk  mbtood. 

flndi*lthOU8h  Ph,M  OM  Pr0,rld*8  *  tha0Patl“^  'bol-proof*  my  of 
finding  a  solution  to  (3)  and  . 

exists  which  i  ’  m±”g  °",t  no  80lutl°n 

It  lntrod  "  4  8en"e*  18  8  solutlon  to  the  whole  problem— still 

it  Introduces  Acuities  of  it.  own.  m  the  first  place,  .xpen- 

certataBt"hOWn  ^  **  aPt  t0  b*  ®dUlT  l0ng  wh,n  p  *“« 

ertain  types  of  structures  mmnmni .n 

ctures,  especially  when  Q  contain,  many  zeros 

tlH  d  "  Wht°h  arl“‘  falrly  °ft#n-  THla  la  ~«T  computing 
Also  introduce,  considerable  round-off  error  Into  th. 

zrms?  M  A°£CL 
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computation.  before  th.  actual  maximisation  of  xQ  la  begun. 

Secondly,  there  la  no  guarantee  that  the  feaalble  baala  obtained  at 
the  end  of  pha.e  one  may  not  be  about  aa  far  from  optimal  aa  poa.lbl. 
The  algorithm  to  be  herein  preaented  la  the  aeoond  of  two  propo.al. 
[8]  for  overcoming  the  aecond  of  theae  objection..  The  flrat  objec¬ 
tion  la  closely  related  to  the  difficulty  mentioned  In  footnote  l— 
that  la,  although  non-convergence  appeara  to  be  extremely  rare,  .low 
convergence  la  fairly  common  for  phaae  one.  However,  It  la  believed 
that  even  thla  difficulty  will  be  Improved  by  the  pre.ent  algorithm. 


111“*  LOOK  AT  ARTIFICIAL  variables  and  the  product  form  of  b-1. 

In  the  uaual  phaae  one,  an  "artificial"  Identity  matrix  1. 
adjoined  to  P  and  artificial  varlablea  are  aaaoclated  with  theae 
unit  vectors,  an  auxiliary  maximising  form  then  being  Introduced 
for  the  purpose  of  eliminating  the  artificial  vectors  or  at  least 
of  driving  the  sum  of  the  artificial  variables  to  sero.  (Thla  aux¬ 
iliary  form  consists  of  a  redundant  equation  Incorporating  an  additional 
variable  similar  to  x0.)  Alao  when  the  product  fora  of  B”1  i.  used, 
thla  la  usually  considered  aa  merely  a  computationally  convenient 
way  of  recording  E T1.  We  now  wish  to  look  at  these  variables  and 
the  elementary  matrices  whose  product  la  E T1  from  .lightly— but 
Importantly — different  viewpoints. 

It  1.  mandatory  that  we  maintain— at  all  stages  of  the  simplex 
process— a  basis  for  (m+l)-apace,  regardless  of  whether  it  1.  wholly 
or  partly  artificial,  or  completely  contained  In  P.  Thu.  when  we 
speak  of  a  "basic,  feasible  »olution"  to  (5)  and  (4): 
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(»)  "solution"  mean.  that  «quatlon  (4)  i.  satlsn.d,  and 
(b)  "feasible"  means  that  condition  (3)  la  satisfied,  whereas 
(0)  "basic"  naans  that  for  all  Xj  >  0.  the  associated  are 
linearly  Independent  (and  Independent  of  Pq)  „d  h.no.  lt 

there  are  p  <  m  such  f_,  (j  >  0)  then  m-p  unit  vectors 
must  "fill  out"  the  basis. 

Clearly  these  three  notions  are  distinct  and  the  Idea  or  a  solution 

being  "basic"  1.  really  ejctran.ou.  to  the  linear  progressing  problem. 
Furthermore , 

(d)  "optimality"  refers  to  (l) 

and  la  still  a  fourth  dlstlnot  notion.  Although  the  simplex  criterion 
for  optimality  requires.  In  praotlee,  the  Inverse  of  a  basis,  the 
degree  of  artlf totality  la  not  a  consideration.  The  essenoe  of  the 
present  Proposal  1.  to  keep  these  four  notion,  distinct  but  to  work 
®n  *11  of  them  simultaneously. 

Since  »  must  always  provide  a  basis,  it  Is  natural  (and  very 
convenient)  to  .tart  with  the  simplest  of  all  bass.,  the  Identity 
matrix  I.  He  will  not  think  of  this  as  an  "artificial  basis"  for 
the  problem  but  simply  a.  part  of  the  necessary  mechanism  on  which 
the  operation  of  the  method  depends.  He  cm  obtain  feasibility 

initially  by  setting  all  Xj  -  0,  j  .  0,1 . .  However,  this 

doe.  not  provide  a  solution  to  (4/  .0  It  is  necessary  to  modify  (4) 

M  that  the  right-hand  aide  1.  null.  For  this  purpose  we  Introduce 
■  auxiliary  non-negative  variables  ^  *  0,  1  -  i,2 . „.8  ^ 

‘It  r  contains  some  unit  vectors  V1  other  than  r,  then  the  x 

for  thtst  vectors  can  be  set  ecru&l  to  h  a  a.  ®  i 

Introduced  for  this  set.  i  ^^lally  and  no  u^  need  be 
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modify  (2)  to 


(6) 


n 


jii  au  xj  "  bi  -  ui 


^  ™  1*2,  a , ,  jin. 


rr: ~ttins  ^ b*  p~vid-  *  ■«-«-  *. 

(°  **.,_*«»„  u.{o,  Ul.  u, . 

“*  then  seek  both  to  maximise  x0  „d  t0  reduce  „  to 

«hn.  "*thtalnlng  feasibility  and  continuing  to  provide  .  bssi.  with 
which  to  operate.^ 

Clesrly  if  P  contains  any  non-null  vector  beside.  t0,  sow 
t}  (3  >  0)  can  be  Introduced  into  the  basi.  ellinatlng  of  the 

"It  '"aiora  V  *  >  0.  «»  choice  of  which  Pj  to  Introduce  win 
be  mad.  to  depend  on  the  optimising  form  (l)  ..  lon8  „  lt  lt  pot_ 

Slbl.  to  improve  ,#.  where.,  the  choice  of  which  0,  to  drop  will  be 

such  a.  to  maintain  feasibility  (l.e.,  keep  (,)  eatlafled)  end  al.„ 

maintain  all  £  0.  The  first  substitution  must  reduce  soma  u 

to  aero  (uni...  certain  u,  were  not  Introduced,  ...  footnote  8) 'when 

the  corresponding  U,  1.  eliminated  from  B.  Whenever  this  haPP«,., 

the  u,  1.  then  not  allows  to  renter  the  problem,  that  is,  when 

“l  ^  *11"lnat*<1  fro“  B«  U1  *•  rsmoved  from  V  and  neither  are  ever 
considered  again.  Ho»ver,  some  u,  may  forwily  remain  i„  0  lt 

lsv.1.  Whenever.  In  the  sequel,  w.  refer  to  2  uk,  the  sum  1.  con¬ 
sidered  as  taken  over  those  uk  still  formally  in  o. 

propoaed^bjMS^ .  ^Bsa le^ and^eiabo  rated #b®*nf  t?tf  byr*  similar  «•*!«. 
varlabla  measuring  tha  sua  of  tha  ^0  1x1  "h*0**  a 

tha  right  mambar  of  tha  radundant  aouatloS  lntroduoad  in 

one.  Nota  that  wa  hava  not^^S^qSISonT  •»»>»»/•<«  in  phase 


i 
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On.  of  the  U±  1.  not  n.o...arlly  eliminated  on  ..oh  iteration; 

“*  oertaln  PJ  b.  Introduced  Into  B  and  replaced  i.t.r 

by  oth.r  Pr  However,  th.  choice  of  a  f ^  to  Introduce  depend,  on 
(1)  and  nof  on  an  auxiliary  for*. 

Baoh  Iteration  i  win  produce  an  elementary  matrix  E.  which 

tran.form.  th.  Inver.,  of  th.  (i-a)*  *  bael.  *  to  th.  Inver.,  of 
the  ith  baele: 

BJ  -  -  BjE^  . . .  *1E<)  ,  H0  -  I  -  ET1  . 

If  ™dwh.n  a  f...ibi.  aolutlon  to  (3)  and  (4)  1.  obtained,  «y  at 
the  N  Iteration,  then  the  ...oolated  bael.  ha.  an  Inver.. 

^  "  *B*H-1  •••  •••  Bj 


where  V*  -  Q  and  for  any  (not  in  ?)  remaining  In  B^  v.  -  0. 
Heno.  <*n  think  of  th.  N  Iteration,  a.  ...entl.lly  computing  the 
Bi  whoa.  product  give.  th.  invere.  of  a  f.a.lbl.  ba.i.  V  Hot.  that 
.  feasible  bael.  may  be  partly  "artificial . "  Cm  the  other  hand  a 
b..l.  choaen  completely  from  00lUBn.  of  p  „ot  provld.  .  ^ 

aolutlon  to  (5), (4).  Hence  th.  notion  of  a  "fea.ibl.  bael."  1.  r,.la 
*  composite  idea. 


XV— REQUIRED  MODIFICATION  OF  THE  F1R3T  SIMPLEX  CRITERION 

Qiolca  A:  (Some  8  <  Q) 

"  » 

There  are  two  main  criteria  which  together  with  Oauaelan 
elimination  ...ontlally  oon.tltute  th.  .Implex  method.  Fir»t,  given 
a  bael.  B,  th.  fir.t  row  of  IT1,  P<>,  i.  ua.d  a.  a  "pricing  vector." 


The  inner  products 
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(8)  ,  J  -  1,2, ... ,n 

are  formed.  If  any  Oj  <  0,  then  an  index  s  is  chosen  by 

(9)  0.  -  min  6.  <  0  , 

8  J  J 

the  usual  convention  being  that  the  smallest  such  index  is  ohosen  in 
case  of  ties.  If  (9)  applies,  then  the  vector  P  is  Introduced  in 

S 

the  amount  O  >  0.  The  new  value  x*  of  the  new  solution,  obtained 
when  P8  replaces  (if  possible)  some  vector  in  B,  is  given  by 

(10)  x*  *  xQ—  66  #  £  xQ  . 

This  same  rule  will  be  used  in  the  present  method  whenever  (9) 
applies,  l.e.,  whenever  some  8^  <  0.  We  will  call  this  selection 
criterion  "Choice  A." 

Choice  Bi  (All  Bj  £  all  u^  ■  0) 

In  the  usual  simplex  procedure,  the  condition  that  all  8^  >  O 
implies  that  an  optimal  solution  has  been  obtained,  namely  the  basis 
B  and  its  associated  solution  vector  V  -  B“1Q.  This  Implication 
rests  on  the  following  two  assumptions; 

(I)  V  has  been  maintained  feasible  on  every  iteration; 

(II)  either  B  contains  no  artificial  vectors  or,  if  it 
does,  the  sum  of  the  corresponding  artificial 
variables  is  zero. 

While  we  will  maintain  feasibility  in  the  extended  sense  that  (2) 
holds  and  also  all  u1  £  0,  we  will  certainly  have  artificial  unit 
vectors  in  B  for  at  least  the  first  m-u  iterations  (where  u  is  the 
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“rc*  “  v  j  > « ~  -  — — 

■.  *  • — -  ~  r  -  r  rrr — - 

i«  equivalent  to  the  aon(H„  ,  *  1  eventuality 

n®  condition  of  optlmalitv 

Prooadura,  and  hanea  ua  y  In  tha  normal  almpla* 

*"°®  ”  tsmlnat.  th.  prooadura  m  anoh  . 

“  °Ptl,,U"'  -^on  ha.  b..„  attalnad.  °“# 

—  Ci  -  'iia-ij  >  o,  so«'_u1-2_o, 

*n  411  °J  ^  w*  °*”  Obtain  no  Improvaoant  in 

™.tmg  faaclbllltjr  (a,.  (l0)J.  but  „  .l<o  “  ;°  ~h- 

*•*  ««»-Pt  to  raduc.  thaw  u  to  ,  *  1  >  «.  than  w. 

solution  to  (4)  p.  1  In  ord#r  to  olalm  a 

r  this  purpose ,  we  borrow  the  »  u  -  / 

“odlfi.d)  from  tb.  dual  .i^lwt  algorlthM_  *  0rlt‘rl0"  <««*«* 

Since  in  fch-i- _ _ _  *  J 


-w — -  L; 

olnoe  In  this  case  we  wi.w  «. 

wl»h  to  reduce  the 

tel  a.  _  _  .  _ 


-  — trm  auce  the  u  so  . «. 
loglo.1  to  sorb  first  on  th.  Ian,,..,  on„  _  1  '  “ 

t  by  *•  Hdno*  «  choose  the  index 

(ll)  ut  -  max  u±  >  o  , 

(taking  the  smallest  sueh  index  say  in 

Ot  to  form  th.  innar  P«du.t.  P  ,  ..7  ~  “* 

; — -  v-—  if  data^rrr^ 

*lth  '•  d°‘*  vlolata  f.Mlbllity_th.n  p  P  .*  r*Pl*°lng  °t 

the  -pivot  element"  for.  *.*  ,  ^  *  yt  ~  yr  **«•  be 

lament  for  the  elimination.  (Since  Y  -  «-lp 

Th*n  al>°»  »»  will  have  »  .  ut/  *  *■'  yl  *  pip 

•o  we  must  have  y  -  a  p  \  ^  *°d  thA*  mu*t  be  non-negativ« 

lnoraaaa  *  a.  7  t!  ‘  •Ithongh  ..  e„not 

o  ^  14  «  ll“l®  «  Posslbl,.  Slnoa  th. 
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change  in  xQ  la  •p&g  -  Y0  where  yQ  ■  P0P,  •  6#  £  0,  we  choose 
•  by  the  following  trlpla  rula. 

(12a)  If  any  PtPj  >  0,  ohooaa  tha  araallast  Index  a  such  that 
yo  -  PoPa  4  Po*! 

rr  "  T  •  v.  >  °- 


(12b)  If  all  ptPj  -  0,  than  oan  navar  ba  raducad  and  thara 
la  no  solution  to  (4).  (Ut  la  Indapandant  of  all  p^.) 

(12c)  If  all  0tPj  £  o  and  soma  Pt?J  <  0,  than  any  allowable 
change  In  B  at  this  point  will  not  deorease  u^.  In  this 
case  we  seek  to  reduce  tha  sum  £  ^  of  tha  ^  still 
formally  in  U.  Forming  tha  corresponding  sum  p£  •  £  Pk 
wa  compute  for  all  J  >  0.  If  any  ^  >  o,  use 
(12a)  to  determine  s  '(replacing  with  p^).  If  all 

PZPJ  ^  °'  th#n  Z  'he  >  0  lB  **«i«um,  l.e.,  cannot  be 
reduced  without  some  x^  going  negative  and  there  Is  no 
feasible  solution  to  (3), (4).  (This  is  equivalent  to 
the  condition  In  tha  normal  phase  one  whan  it  Is  deter¬ 
mined  that  no  feasible  solution  exists.) 

Tha  reader  will  note  that  (12c)  Is  nothing  but  a  different  way  of 
operating  the  usual  phase  one.  However,  this  rule  does  not  eome 
Into  play  until  all  else  has  failed.  Whatever  novelty  the  present 
proposal  can  claim  lies  In  this  change  of  order  and  emphasis  in 
tackling  the  four  requirements  (a),  (b),  (c),  (d).  Although  the 
rules  of  selection  for  vectors  to  enter  In  and  drop  from  the  basis 
are  more  complicated  here  than  In  the  norsml  procedure,  this  Is  a 
secondary  consideration  compared  with  time  for  and  accuracy  of 


r  » 


~ii- 


•olutlon  on  a  largo  alaotronlc  oomputer. 

If  (12)  allows  a  P  to  ba  ahm.n 
It  in  t.™  of  B  “  U"“1  ” 

(13) 


"  B"1P*  ‘  {'».  r,*  ....  y„} 


7  U  n°  S“r“t**  th*‘  «t  >•  eliminated  «  feasibility 

2TT*  ""  ",OU8h  (l2a)  what^v8r  change  i. 

m  T  *nt*r  B  Wlth  “*  °°,mol*nt  V*iO.  The  o^g. 

in  ut  (or  2  ^  if  (i2c)  wa.  u.#d)  W±1I  b# 

ut  "  ut  -  ^  ut 

*°r  2  “k  "Iuk-#Jyk^2uk  In  case  (l2c))  . 

Again  It  should  b.  emphaalaed  that  whan  Choice  C  auet  be  used 
oon.id.«bly,  then.  1.  no  guarantee  that  the  pre..„t  Mtho<J  wU1 
y  fa* ter  th«  the  non»l  phae.  one.  However,  1,  then.  1.  .  oh.no. 

~~~  “  ^ ione  m-b*r  °f  «««~.  p«~»t 

«thod  leave,  the  door  opan,  a.  it  a.,,.  for  auoh  .  abort  out  to  be 


are 


_  0  “  11  P*fh,p*  worth  “•"tlonlng  that  when  moat  of  the  a  . 

ft.««lr  *1‘llk*ly  t0  b*  »—  the  outset  and  for 

Ions.  On  the  other  hand.  If  most  (or  all)  th.  .  y  •  ‘ 

P~onsCa^tMr^ol?eJ*JSrIt.ff;rSU^*tiSd  P#rhap8  tile 

thT^Jf  °f  —  S^nW^-ob^iHf 
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V— MODIFICATIONS  OF  THE  SECOND  SIMPLEX  CRITERION 

Th«  ••cond  criterion  In  the  simplex  method  la  oonoarnad  with 
tarmlnlng  the  Index  r  of  the  baala  veotor  to  be  replaoed  by  P 
Representing  Pg  In  terms  of  B  by  (13),  the  usual  criterion  la  to 
form  the  ratios  (provided  any  y^  >  0) 

vi 

(14)  ~  for  all  y±  >  0  and  1  >  J 
and  then  to  choose  r  by 

(15)  *r  -  min  for  all  l  lnoluded  In  (l4), 

whence  *r  >  0.  Providing  not  all  y±  £  0  this  same  orlterlon  will 
be  vised  here. 

In  the  usual  procedure,  the  case  of  all  y±  ^  0  (l  >  0) 
that  xG  has  no  finite  maximum.  For  the  present  method, 
however,  this  last  lnq>lloatlon  Is  no  longer  valid.  It  must  be  modi¬ 
fied  and  for  this  purpose  we  distinguish  three  eases.  We  first  note 

that  all  yA  <;  0  can  only  oocur  in  the  oase  0^  <  0.  (See  Choloe 
A  In  IV.) 

Case  I.  All  u^  ■  0,  all  y ^  £  0,  and  yk  ■  0  for  all 

remaining  In  U.  In  this  oade  a  olass  of  feasible  solutions 
*o  (5), (4)  oan  be  oonstruoted,  whose  values  have  no  upper 
bound,  namely, 

(16)  B(V-dYs)+gPB 

with  a  new  value  x*  given  by  (10).  Slnoe  y„  -  a  p.  -  »  ‘  <  o 

w  '  O  rO  8  8  * 

dearly 
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«»••  equivalent  to  the  o*$t  In  the  normal  simplex 
mathod  (phase  two)  where  all  y±  £  0  and  tha  iterative 
prooadura  tormina taa. 

Hie  function  of  tha  redundant  aquation  during  phaee  two  of  tha  nonaal 
prooadura  la  to  inaura  that  an  artificial  variable  does  not  gat  at?va 
zero  alnca  this  would  require  tha  redundant  variable  to  go  below 
zero  (or  vice  veraa)  which  la  prohibited  by  (14)  and  (15).  Since  in 
the  preaent  method  no  redundant  aquation  la  used,  thia  poaalblllty 
inuat  be  avoided  by  examining  the  yk  for  whioh  la  in  U.  Thia  leada 
us  to  oaae  two. 

£•••  II*  All  *±  m  0,  all  y^  ^  0,  but  aome  yk  <  0  for 

\  in  u*  Hare  Pg  can  be  Introduced  into  B  by  eliminating 
one  of  the  in  B  but*  not  in  P.  For,  since  all  u^  •  0, 

*  Ur/yr  •  0  (yp  <  0)  and  the  new  aolutlon  vector  V* 
will  not  differ  from  tha  preaent  aolutlon  V  (of.  (l6) )  and 
hence  faaalblllty  la  maintained  whJ la  one  more  u^  la 
dlaearded. 

Case  III.  All  y±  £  0  but  aome  u^  >  0.  (Again  note  that 
»,  <  0.)  If  >  0,  than  x0  will  lnoreaee  but  ao  will 
all  ^  for  whioh  yk  <  0.  Since  there  la  no  point  in 
inoreaalng  xQ  indefinitely  at  the  expenee  of  llkewiee 
increasing  deviatlona  from  Q  on  the  right-hand  aide,  we 
arbitrarily  reject  the  chosen  ?a  at  thia  point  and  return 
to  Choice  C  in  eeotlon  IV.  In  thia  case.  Choice  C  may 
actually  result  in  an  increase  in  xQt  but  for  the  present 
iteration  we  are  concerned  only  with  reducing  £  u^  of  the 
still  formally  in  U. 


RM--1275 


RSPKRKNGBS 


[l]  Dantsig,  George  B.,  Alex  Orden,  end  Philip  Wolfe,  "the 
Oenerelised  Simplex  Method  for  Minimising  M  Id 
Form  under  Linear  Inequality  Restraints,  1 §|ffl 
5  April  1954  (Revised ) . 

[2I  .  .  . . . ,  "Duality  Theorems, *  /Kfi 

30  October  1953  (Revised). 


Method,"  RAID] 


itional  Algorithm  of  the  Simplex 
56,  26  October  1953  (Revised). 


3  ■  Algebraic  Proof  of  the  MinMax  Theorem," 

RARo  Rjt42S£  vUr vised)  18  December  1953. 

-  -  —  .'7  and  Win.  Orchard— Hays ,  "Alternate  Algorithm 

for  the  Revised  Simplex  Method  using  a  product  form 
for  the  inverse,"  RAMDlR^Kj^  19  November  1953. 

. .  -t,  ■  . . '• . -  ' — . —  ,  "The  RAND  Code,  for 

the  Simplex  Method  (Preliminary  information  on  Ooda^ 
OHS  XI  for  the  IBM  701  Electronic  Computer),"  RAMP /RM-12< 
20  January  195*.  I 


- ,  "the  Dual  Simplex  Algori thm, "  RAMP  Rfc-1270 

3  May  1954  (Revised).  ""V].., 

.  . .  "Composite  Simplex-Dual  Simplex  Algorithm — I," 

RAMDRIUm  26  April  1954. 


